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Abstract 

We consider a 2D rotating Bose gas described by the Gross-Pitaevskii (GP) theory and investigate 
the properties of the ground state of the theory for rotational speeds close to the critical speed for 
vortex nucleation. While one could expect that the vortex distribution should be homogeneous within 
the condensate we prove by means of an asymptotic analysis in the strongly interacting (Thomas- 
Fermi) regime that it is not. More precisely we rigorously derive a formula due to Sheehy and 
Radzihovsky [Phys. Rev. A 70, 063620(R) (2004)] for the vortex distribution, a consequence of which 
is that the vortex distribution is strongly inhomogeneous close to the critical speed and gradually 
homogenizes when the rotation speed is increased. 

From the mathematical point of view, a novelty of our approach is that we do not use any 
compactness argument in the proof, but instead provide explicit estimates on the difference between 
the vorticity measure of the GP ground state and the minimizer of a certain renormalized energy 
functional. 

MSC: 35Q55,47J30,76M23. PACS: 03.75.Hh, 47.32.-y, 47.37.+q. 

Keywords: Bose-Einstein condensates, quantized vortices, Gross-Pitaevskii energy, renormalized 
energy, vortex lattices. 

Contents 

1 Introduction 

2 Statement of the Main Results 

2.1 Sketch of Proofs 

2.2 Discussion and Extensions 

2.3 Plan of the Paper 

3 The Renormalized Energy 

4 Energy Upper Bound 

4.1 The Trial Function 

4.2 Evaluation of the reduced energy £ [ftriai] 



1 



CoRREGGI, Rougerie - Inhomogcncous Vortex Patterns in BECs 



2 



5 Energy Lower Bound and Convergence of the Vorticity Measure 

5.1 Preliminary Steps 

5.2 Evaluating the Individual Energy of Vortices 

5.3 Evaluating the Interaction Energy of Vortices 

5.4 Completion of the Lower Bound Proof 

5.5 Convergence of the Vorticity Measure 



mi 

M 
[3D 

[33] 



Appendix 



[37] 



1 Introduction 

Vortex nucleation in equilibrium states of rotating fluids is a hallmark of superfluidity and its experimental 
observation in rotating Bose gases was a key step in the understanding of the properties of Bose-Einstein 
Condensates (BECs). Among other spectacular observations, that of triangular lattice^] containing up to 
hundreds of vortices (see, e.g., [BSSDl lARVKl IMCWDl IRAVXKl ICHESQ gave a new strong motivation 
for theoretical studies. 

A rather natural question, that has not been addressed immediately after the observation of vortex 
lattices, can be formulated as follows. The BECs that are produced in laboratories are trapped gases, 
meaning that the confinement against centrifugal forces is provided by a magneto-optical potential. As 
a consequence the matter density profile of the condensate is not homogeneous and depends strongly on 
the type of trap that is being used. In this respect, it is rather striking that the vortex lattices observed 
in experiments seem to be perfectly homogeneous with a uniform mean distribution of vortices in the 
samples. Indeed, the vortex density could depend on the underlying matter density in various ways, for 
example regions of low matter density such as the boundary of the fluid could attract and pin the vortices. 
How come that such effects do not seem to be observed and that the vortex lattices are homogeneous, at 
least to a very good approximation? 

To our knowledge, this question has been formulated and answered first by Sheehy and Radzihovsky 
in [SR11 ISR2j . More precisely, the relation between the matter density and the vortex density has been 
elucidated based on formal arguments (see also [BPGWJ), leading to a formula whose efficiency has 
been favorably compared to experimental data CHES and to numerical simulations [Dan]. Among the 
findings of |SR1[ ISR2] is the fact that the vortex density does depend on the matter density, but in a 
subtle way that has leading order effect only close to the critical speed for vortex nucleation. It is thus 
not surprising that the vortex lattices seem completely homogeneous in experimental situations since, 
in order to observe many vortices, the rotational velocity is taken well above the first critical speed. 
However, some slight inhomogeneity of the vortex lattice survives for these large angular velocities, as 
predicted by [SR1[ ISR2j , and it can in fact be observed as a small correction to almost uniform vortex 
densities |CHESl iDan] . 

Although it is experimentally difficult to observe the transition regime where the vortex lattice is 
expected to be inhomogeneous, the theoretical question remains of interest. The main model used for 
the description of rotating BECs is the so-called Gross-Pitaevskii theory, which can be rigorously derived 
from the underlying many-body problem (see [LSSY1 ILS] and references therein) in a suitable limit. It 
is of importance to be able to rigorously derive a formula for the vortex density from CP theory. 

In this context, the critical speed for vortex nucleation has been rigorously computed in |IM1} lAJRj 
and the distribution of the first few vortices to appear in the condensate studied in [IM2J. On the other 
hand, the regime well above the first critical speed has been treated in |CY| ICPRY1| ICPRY3) where it 
has been shown that the vortex density is homogeneous to leading order, in the sense that many vortices 
are packed in the condensate and their average distribution is uniform. Note that the latter contributions 

ll Abrikosov lattices' of vortices analogous to those occurring in type-II superconductors. 
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concern the average distribution of vortices in a regime where they are densely packed in the fluid, it does 
not give access to the precise pattern formed by the vortices. A rigorous proof starting from GP theory 
that the vortices arrange on a triangular lattice seems to be still out of reach, despite recent advances in 
the related Ginzburg-Landau (GL) theory |SS4j . 

Summing up, there is still a gap in the rigorous theory of BECs between |IM11 IIM21 [AJR] and 
|CY1 ICPRYll ICPRY3] . This gap corresponds to the regime where the rotation speed is larger than 
the first critical speed but of the same order of magnitude and this is precisely the regime where the 
inhomogeneity of the (average) vortex distribution should come into play The present paper aims at 
filling this gap by rigorously deriving from GP theory the formula for the vortex distribution of SRL fSR2] . 

Our mathematical setting is the following: we consider a two-dimensional rotating BEC confined by a 
trapping potential V(r) = r s (wither = |r|), in the framework of the GP theory. After a suitable scaling 
of length units (see |CPRY31 Section 1.1]), the GP energy functional can be written 

£ GP [y]:=J dr ji|V#| 2 -Q**L* + ^|*| 2 + ffl-j, (1.1) 
where f2 is the angular velocity, 

V(r):=r s , s > 2, (1.2) 

L stands for the third component of the angular momentum, i.e., in polar coordinates r = (r, L = —idn 
or equivalently L — r • V^, V := (—d y , d x ), and the coupling parameter e > is going to be assumed 
small (e -C 1), i.e., we study the so called Thomas-Fermi (TF) limit of strong interactions. The wave 
function "J : M 2 H> C belongs to the domain 

^ GP := G H\R 2 ) n L 4 (M 2 ) : r s |lf G L^R 2 ), H*^^) = l} (1.3) 

and the ground state energy of the system is obtained by the minimization of £ . 

E GP := inf £ GP m. (1.4) 

We denote by <]/ GP any associated minimizer (there is no uniqueness in general). 

We will mostly be interested in spotting the vortices of fy GP , that is its zeros that carry a non-zero 
phase circulation (i.e., non-trivial topological degree or winding number). In particular, connecting to 
the preceding discussion, we would like to derive a relation between the distribution of vortices and the 
matter density of the system. The latter, given by l^ 013 ! 2 , can be approximated by minimizing the 
simplified TF functional 

S TF lp]:=e- 2 [ dr [r s + p) p, (1.5) 



obtained by dropping the kinetic terms in (1.1). Here p > plays the role of the matter density, 



normalized such that J R2 p = 1. The minimizer of (1.5) is the explicit radial function 



p TF (r)= l [A TF_ rS]+; (16) 

where [ • ]+ stands for the positive part and A TF is a normalization parameter that ensures ||p TF || 1 = 1- 
Note that p TF has compact support in the ball of radius R TF and a rather simple computation yields 



ffTF- , X TFf(o+V/> flTF/ A TF)Vs ^F = / 2(s + 2) | 

4(s -I- l)e 2 v ; K J > \ - 1 / 



s/(s+2) 



2 Throughout the paper we use the following convention: a vector will be denoted in bold fonts (e.g., r), whereas normal 
fonts will always denote scalars (e.g., r = |r|), which might also be the modulus of a vector. 
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with E TF standing for the TF ground state energy. Essentially, T> TF is the region occupied by the 
condensate: we are able to prove that the mass contained in E 2 \ T) TF is extremely small, the reason 
being that \I> GP decays exponentially (both as a function of r and e) in this region. 

To spot the vortices of ^> in the bulk T> TF of the condensate, it is standard to consider the so-called 
vorticity measure \i associated with 5' GP . As we are dealing with a regime where the condensate contains 
a large (actually oc | loge|) number of vortices, fi should be interpreted as giving the mean distribution 



of vortices. With the definition we will adopt below (see ( 2.13 1), one can actually prove that, if W 



rGP 



contains J vortices of degrees d±,...,dj and locations ai , . . . , a,/ 

J 

H^2K\\oge\- 1 Y J d J 5 BLj (1.8) 

i=i 

in the (C*)* topology, where 8 &j stands for the Dirac mass at a.,-. Note the scaling factor |loge| _1 in 
(1.8): since the condensate contains (D(\ loge|) vortices, n will be a quantity of order 1. 

The mechanism for vortex nucleation in rotating superfluids is now well understood, see, e.g., t AAB. 
IIM1) IIM2| ICRYj . A vortex becomes favorable in the system if it can lower the energy by interacting with 
the rotation field. More precisely, the interaction with the rotation field should overcome the energetic 
cost for vortex nucleation, given by 

7r|log £ ||^|p TF (a i ) (1.9) 

with dj the degree of the vortex and &j its location. To evaluate the energy gain brought by the vortex 
compensating the rotation, it is convenient to introduce the potential function 

F TF (r) = --^ [ RT \ttp™(t). (1.10) 
I logej J r 

The energetic gain of a vortex is then 

2 7 r|loge|d JJ F 1 (a J ). (1.11) 

By comparing the cost and gain of vortex nucleation it is obvious that in order to be energetically 
favorable, vortices should have positive degrees dj, because F TF is clearly negative. These considerations 
lead to the definition of a radial function giving the energetic cost of a vortex of degree 1 located at &j 
with |aj| = r: 

ffTF (r)= l p TF (r) + F TF (r) . (1.12) 

Looking for its minimum, one finds that it lies at r = 0, indicating that this is where a vortex is most 
favorable. Equating gain and cost of a vortex at the origin, one can see that the critical speed for vortex 
nucleation is given by 

nc^nxlbgei, n i: =£(f^±4) , (Lis) 

2 \ ITS J 

namely _ff TF (0) > for il < fl ci and H TF (0) < for > f2 ci . Note that the above value is that one 
finds by applying the rigorous analysis of jIMll lAJRj . 

The question we address in this paper is what happens when Q is chosen of the form 

n = n |io g e|, n >fii, n = o(i), (1.14) 

i.e., when f2 is strictly larger than the first critical speed but of the same order of magnitude when e — > 0. 
We prove that in this regime the vorticity measure fi is to leading ordeij^] 



l{ffTF< 0} (1.15) 



3 We denote by I5 the characteristic function of the set <S. 
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in a sense to be made precise in the next section. Computing the term Vi? TF this can be rewritten 

as 

/i(r)« i^log(p TF (r)) +2f> 1 {h tf< 0} , 



(1.16) 



which is the formula found by Sheehy and Radzihovsky, except that our analysis shows that vortices 
should lie only where the cost function is negative, which was not clearly mentioned in |SR2| . Note that 
vortices are also confined to the region where 

1 -d 2 r \og(p TF (r))+2n >0, 

which is far from obvious if only the cost and gain considerations sketched above are taken into account. 



Indeed, the particular form (1.16 1 is to a large extent due to the interaction between vortices. 



Several interesting properties of the (average) vortex distribution as a function of Q can be read off 



from ( 1.16 1 



When fir, < fii in (1.13), a straightforward computation reveals that the cost function H is 



positive everywhere. Vortices are thus not favorable and the vortex distribution vanishes identically. 
We thus recover the expression of |IM1[ lAJRj for the first critical speed. 



In the regime (1.14), one can compute (see Section [3]) that there is a non-empty region where 
H TF < and V^4f Vi? TF > 0, whose size increases with increasing fio until it finally fills the 
whole sample in the limit fio — > oo. This region is, according to (1.15), filled with vortices, a 



behavior that is reminiscent of the 'obstacle problem regime' in GL theory [SSI . 



The first term in the right-hand side of (1.16) can not be constant, except when p TF is constant 
itself, which can only happen in the somewhat unrealistic case of the flat trap considered, e.g., in 
[CPRY1 . This term is thus responsible for an inhomogeneity of the vortex distribution, whereas the 
second term yields a constant contribution of 2f2o (20 in the physical variables) units of vorticity 
per unit area. 



• The first term in (1.16), responsible for the inhomogeneity, is independent of Oo- Its importance 
relative to the second one thus diminishes with increasing Oo. The inhomogeneity then becomes 
a second order correction in the limit O — > oo, which corresponds to O S> |loge|. As discussed 
above, this is observed in experiments and numerical simulations. This also bridges with the 
situation considered in [CPRY3] . where we proved that, if O ^> | log e| , the vortex density is to 
leading order constant and proportional to 20. 

The next section is devoted to a more precise statement and discussion of our results. Proofs are 
jiven in Sections [3j [4] and [5] 
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2 Statement of the Main Results 

We now turn to a more precise description of our results. Recall that we consider the asymptotic behavior 



of the ground state energy and minimizer of the functional ( 1.1 ) when e — > with the scaling 

fi = fi |loge|, n >f2i, 



(2.1) 



where f2o is assumed to be constant for the sake of simplicity and Qi is defined in (1.13). The regime 



^o^l corresponds to what we have considered in |CPRY3j , while in |IM11 IIM2 it was studied the case 
Q = f2 1 |log£| + 0(log|loge|). 

To take into account the inhomogeneous density profile it is convenient to introduce the following 
energy functional 

£ GP [/] 



J dr |I|V/| 2 +e- 2 [r* + / 2 ] / 2 | 



with ground state energy 



£ UP := min £^[f] 

/e@ GP 



/ = /*} 



(2.2) 



(2.3) 



Standard arguments show that there is a unique strictly positive radial minimizer which will be denoted 
by g. Note that £ GP coincides with £ GP restricted to real functions. In this sense, its minimization 
corresponds to the search for a vortex-free profile. The main difference between £ GP and £ TF is that the 
former includes the contribution of the radial kinetic energy due to the bending of the density profile. 
We will prove in the Appendix that g 2 « p TF in a suitable sense. 



The expression (1.15) enters our problem through the minimization of a 'renormalized energy' (we 



employ the consecrated terminology of GL theory BBH, SS2J) expressing the energy of a given vorticity 
measure v in terms of the TF density p TF in units of | logej 2 



X"[„] 



1 



„TF 



|V/l„| 



1 



p tf W\ 



wher^f] 

and 



TF •=supp(p TF )=6(i? TF ), 



F TF v 



V 



-V(4fV/i„) =finI' TF , 



on dV TF . 



(2.4) 

(2.5) 
(2.6) 



Recall that F TF is defined in (1.10). The minimization of the renormalized energy in its natural energy 
space 



t> tf L P 



1 



\VK\ 2 + P 



TF I 



< +oo 



(2.7) 



is discussed in details in Section pi We prove (see Theorem 3.1) that X TF has a unique minimizer p* 
among the measures in ^# p tf (V tf ) . It is explicitly given by 



p-k 



n TF 



VH 



TF 



1 



{H TF <0}; 



(2.8) 



4 In all the paper, B(q) denotes the disc centered at the origin with radius g, while B(r, q) is the same disc but centered 
at r £ IR 2 . Also we will sometimes omit the measure in the integral (as, e.g., in the first term of l |2.4[ l), when it is the usual 
two-dimensional Lebesgue measure dr. 
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i.e., it is exactly the measure appearing in the right-hand side of (1.151. We also set 

1 



7-TF 



I'* 



supp(/J* 



where the second equality is proved in Section [3] below. Note that by (2.8) one easily has 



/ TF < 0, 



(2.9) 



(2.10) 



since /i* > and H TF < on the support of /i*. Moreover both the renormalized energy as well as its 
minimizer /i* are fixed when e — > 0, thanks to the extraction of a scaling factor | loge| 2 . 
We can now formulate our first result about the GP ground state energy asymptotics: 



Theorem 2.1 (Ground state energy asymptotics). 

= Oq | log e| , with Hq > Vti, then 



E 



GP 



E 



GP 



7 TF |log^ I 1 + 



log I logs I 
|loge|V2 



(2.11) 



in the limit e — > 0. 



Remark 2.1 (Composition of the ground state energy) 

The leading order contribution to E GP is given by E GP , which in turn contains E TF of order e~ 2 (see 
(1.7)) and a remainder due to the radial kinetic energy of the vortex- free profile. Such a correction can 
be shown to be of order | log er | CPRY3, Proposition 2.1], i.e., much smaller than the contribution of 
vortices I TF | loge| 2 , which on the other hand is a rather small correction to the main term E TF . 

□ 

Remark 2.2 (The renormalized energy) 



Note that the functional (2.4) is defined for the 'reasonable' vorticity measures in ^ p tw(V ), in par- 
ticular those arising from wave functions in the manner of (2.13) below. It is not well-defined for sums 
of Dirac masses, which is a significant difficulty in the analysis, in particular in view of (1.8). The first 



term in (2.4) corresponds to the interaction between vortices, it is computed in a way reminiscent of 
electrostatics: h v is similar to a potential generated by individual electric charges distributed according 



to the charge density v and the first term in (2.4) is the corresponding electrostatic energy. The other two 



terms can be understood as the sum of the cost and gain due to each individual vortex, in view of (1.8), 
(1.9) and ( |1.11[ ). Note that, if it was a priori known that the minimization could be restricted to positive 
measures and the density p TF was constant, (2.4) would reduce exactly to the electrostatic energy of a 
positive charge distribution in the potential \p rt + F TF . In this case (2.6) would become the Poisson 
equation for the charge distribution v. h v could then be interpreted as an electrostatic potential and its 
gradient as the corresponding field. In this analogy, the non-constant weight l/p TF can be thought of as 
modeling a sample with non-homogeneous conductivity. 

□ 

As we will prove below, the minimization of £ gives the matter density of the system to a very 



good approximation, i.e., |\& 



GP|2 



g 2 . It is thus natural to write ^P GP in the form 



(2.12) 



where u is essentially a phase factor accounting for the vortices of >I> GP , times a profile vanishing close 
to the vortex cores and almost equal to 1 elsewhere. A convenient way of spotting the vortices contained 
in u is to use the so-called vorticity measure 



[i := | loge| x curl 



(uVu* - m*Vu) 



(2.13) 
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which is (up to the | logej 1 factor) nothing but the curl of the superfluid current 



(uVu* — ti'Vu) 



(2.14) 



and thus (in analogy with fluid mechanics) a good candidate to count the vortices of u. 

As usual [AAB , IM2] |R] , energy methods do not allow to spot vortices lying too close to the boundary 
of the domain. Indeed, p TF vanishes on <92? TF and thus, according to (1.9), a vortex lying close to the 



boundary carries very little energy. We will thus limit ourselves to analyze the behavior of /i in the 
smaller ball S(i?buik) with i?buik satisfying 



-Rbulk < R 



TF 



\R 



bulk 



R TF \=o{n- 1 ). 



(2.15) 



In fact in (5.11) below we will make the precise choice i?buik = R TF — CO, 1 for some given explicit 



constant C . However our proof works just the same provided the fixed constant C is chosen small 
enough. 

Note that by restricting ourselves to i?(i?buik) we are om y neglecting a small part of T> TF and the 
domain B(-Rbulk) contains the bulk of the mass of the condensate, in the sense thatQ 



dr|* GP | 2 = l-o(l) 



(2.16) 



B(flbuik) 



in the limit e — » 0. 

We prove that fi is close to /x+ in S(-Rbulk)> with the meaning of 'close to' specified by the following 
norm, which is defined for measures v. 



sup 

^eCi^-Rbuik)) 



/ 



dr 



1 

Ttf 



1/2 



(2.17) 



|V0| 2 



\JB(fl bulk ) P 

Theorem 2.2 (Asymptotics for the vorticity measure). 



|V0|| 



L~(B(fibuik)) 



Let /i and /i* be defined respectively in (2.13) and (2.8). Then we have 

'log | logel 1 / 2 " 



11^-^*11 



< O 



l0g£ 



1/4 



(2.18) 



in the limit e — > 0. 



Remark 2.3 (The norm || • || TF ) 



The expression (2.17) is a rather natural definition for the norm of a measure. The norm of the test 
function (j> appearing in the denominator contains two parts. The first is naturally associated with the 



energy functional (2.9), so an optimal statement must necessarily include this term. The second term 
|| V0||ioo(g(^ b \\ appears when regularizing fi in the course of the proof. Without this term, the norm 
||-|| p tf would not be well-defined for a Dirac mass, and in view of (1.8) this would be rather problematic 
in our setting. As we will see in the proof, a suitable regularization of /i can be estimated in the norm 
where the term ||V0||£°o/g( flbiill \) is removed in the denominator of (2.17). We still believe that it is 



necessary to include that term to state a result on the asymptotics of /i. 

Note finally that if p TF was uniformly bounded below by a positive constant independent of e in 
£>(i?buik) (which is not the case), the norm (2.17) would be equivalent to the (C* (B(-Rbulk))) norm. In 



5 This is an easy consequence of the fact that |*E' GP | is uniformly bounded (see Appendix). 
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particular, this means that for any fixed R < R TF and e small enough, (2.18) yields an estimate on the 
{Cl(B(R))) norm of fi— Indeed, for e sm all enough, B(R) C S(-Rbulk) an( i P TF is bounded below in 



We can thus deduce from Theorem 2.2 that for any fixed R < R 



? TF 



pi -> M* 

strongly in (C£(B(.fl)))* as £ ->• 0. The reason why we consider the approximation of /j, by ^+ in the 
larger ball S(-Rbulk) has already been explained: such a ball contains the bulk of the mass in the e — > 
limit, whereas does not. 

□ 

Remark 2.4 (Asymptotics for an explicit vorticity measure) 

It is also possible to state the result in terms of an 'explicit' vorticity measure of the form of the right- 
hand side of (1.8), in the spirit of, e.g., [R] Theorem 1.2] or |CPRY1[ Theorem 1.1]. Indeed, the proof of 
Theorem 2.2 requires to localize the possible vortices in small balls of centers &j and radii Qj, j — 1 . . . J, 



and any statement on the vorticity measure p translates into one on the measure 



2,7 log , ^2<1 



thanks to the so-called Jacobian estimate (see Proposition 5.3 below). Note however that the norm in 



which the estimates hold in this case is necessarily weaker since Dirac masses are less regular than p. 

□ 

We will comment further on these results below. Section |2.1| presents a sketch of our proofs and 
Section |2.2| contains a comparison to earlier results and discusses a possible extension of our method to 
a different setting, namely that of the third critical speed in a flat trap, studied before in (CRY! IR] . 



2.1 Sketch of Proofs 



For the convenience of the reader we now sketch the main ideas of our proofs. Several standard techniques 
will be employed and some of the ideas we use originate in jABM| and [R] (one can in particular compare 
the following sketch to [R] Section 1.2]). We focus on the way the energy lower bound is obtained because 
we believe this is more helpful in explaining the origin of the renormalized energy (2.4). We also slightly 



deviate from the actual proof procedure in several places when this serves the purpose of our heuristic 
considerations. 

As is standard, the proof starts with an energy decoupling: 



E GP = £GP + 



dr 



V|Vu| 2 - g 2 ttv^ ■ (iu, Vu) + 9 
2 e 



2 (i-M 2 ) 2 



(2.19) 



where (iu, Vu) = j is the superfluid current defined in (2.14). The only input in (2.19) is the variational 
equation for g. Then, as already mentioned, the bulk of the mass is contained in 2? TF , so that we make 
a very small error by restricting the integration to this domain. Also, since g 2 f=a p TF , we can, at least at 
the level of heuristics, consider the reduced functional 



S[u] 



drl -p 

■pTF [ 2 



TP 



p^Qr- 1 - ■ (iu,Vu) 



P 



TF Z 



-(i-M 2 ) 2 



(2.20) 



and we essentially have to understand why its minimization reduces to that of I TF , once a scaling factor 
0(| loge| 2 ) has been extracted. It is fairly easy to obtain from the rough upper bound (trial state v = 1) 



£[u] < 0, 
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that the sum of the first and the last term in (2.20) are suitably bounded above. Using the vortex balls 



method introduced independently in [Jel[Sa], this allows to control the area of the set where \u\ differs 
significantly from 1 and enclose it in a finite collection of disjoint small balls Bj, j = 1, . . . , J that serve 
as 'approximate vortices'. One also proves that each ball contains a kinetic energy whose leading order 



is precisely the energetic cost we have been alluding to in (1.9) 



jf dr ip TF |V^| 2 > 7r|^|p TF (a,) |lo gi 



(2.21) 



where dj is the degree of u on dBj and a.j the center of Bj . 

To evaluate the energetic gain of vortices it is convenient to integrate by parts the second term of 



(2.20) using the potential defined in (1.10): 



dr p^r 1 - ■ (iu,Vu) = |loge| / dr F TF curl(w, Vu). 



(2.22) 



Then according to (2.13) and (1.8) (which can be put on a rigorous basis thanks to the Jacobian estimate 
method [JS2] 1 

(2.23) 



loge| / dr F TF curl(m, Vw) « 2tt| lo g£ r| V^F TF ( 



Note that the expressions for the energetic cost and gain of vortices can be motivated by a simple 
computation involving an ansatz of the form 

where we have used complex notation z := x + iy for a point r = (x, y) £ R 2 and Zj are the positions of 
the vortices. The real- valued function £ is a cut-off ensuring that u vanishes close to the vortices. What 
actually comes out of such a computation is a factor <i 2 instead of \dj\ in ( |2.21 1 but rigorous analysis has 
so far been limited to the obtention of the smaller \dj\ factor (recall that dj is an integer), except when 
a priori information is available on the vortex distribution. The difference is anyway of no concern to us 
because it is always favorable for the vortices to be singly-quantized, i.e., dj = 1. 

There now remains to bound from below the part of the kinetic energy contained outside the vortex 
balls. To this end we note that, since outside UjBj we have \u\ 1, 



P TF \Uj8. 



dr^ TF |V M | 2 



X> TF \U.|Bj 



dr^ TF |j| 2 



with j defined in ( 2.14 ). Now comes the procedure of regularization of the vorticity measure we have been 



alluding to in Remark |2|2| since the vortex balls contain vortices, i.e., phase singularities, the gradient 
of the phase of u is expected not to be very well behaved inside them. We thus exclude these regions by 
setting 

f.i, in V TF \ UjBj 



J := 



0, in Uj Bj, 



(2.24) 



and write 



1 



dr ;p TF |Vu| 2 



X> TF \U 3 Bj 



TFi~|2 



dr -p- |J 



With j is associated the vorticity measure (rescaled as p.) 

p := |loge| _1 curlj 
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and in turn we can associate with fl the potential defined as in (2.6). We then have by definition 



1 



curl ( j + | log e\ -^ ¥ W ± h jl ) = 0, 



(2.25) 



which implies that j and | logej^p V 1 ^ differ by the gradient of an H 1 function. Then we conclude (see 



Lemma 5.2 for details) 



/ drV F |V W | 2 >i|log £ | 2 / dr^plVfy 



(2.26) 



and this term has to be interpreted as the energy due to the interaction between vortices via the potential 
hjx they create. 

Now, as the term 'regularization' suggests, one can actually prove that fl /i in a suitable sense akin 



to that in which (1.8) can be made rigorous. Gathering (2.21), (2.22), (2.23) and (2.26) and dropping 



the last term in (2.20) we may thus write (neglecting all remainder terms for simplicity) 



£[i 



> 



\oge\ 



•pTF 



1 



2p 



TF I V M 



i 



p tf IAI 



F TF jl 



(2.27) 



where we have again used the informal relation (1.8 1. This is of course the desired lower bound, since 



one can recognize in the right-hand side our renormalized energy functional. 

Now come the main technical novelties of the paper. Observe first that we do not prove a priori as 
is often done [AABl IABM1 IIM21 IR] that most vortices carry a positive degree. This fact follows directly 
from the minimization of I TF . In other words, the minimizer of I TF is automatically positive and there 
is no need to restrict the minimization to positive measures. Indeed, although the problem of minimizing 
I TF bears some resemblance with the obstacle problem in GL theory |SS21 Chapter 7], it is in fact a bit 
simpler and its explicit unique minimizer is given by (2.8). Contrary to the GL case we do not obtain 



a free boundary problem: the region where vortices should lie is directly encoded in the cost function 
I TF has the following very nice stability property (see Section [3] for the proof): for any 



(1.12). 



Moreover 
measure v 



z TF H > i TF 

We can thus conclude from the above that 



1 



■pTF 



2p 



TF 



\Vh v - 



(2.28) 



£ GP >£ GP +/ TF |loge| 2 + |log e | : 



and the last term is a norm squared of /2, corresponding to 



•pTF 



2p 



but with the ||V^>||l°° term dropped 



in (2.17), as one case easily see by a simple duality argument. 



What remains to be done is the construction of a trial state giving an energy upper bound confirming 
that the above lower bound is optimal. We adapt a well-established technique (see, e.g., |AAB[ [ABM, 
ISS21 |R] ) based on a Green representation of the potentials defined as in (2.6) to out setting. We use a 
Riemann approximation of /i* by a measure located in small balls that mimic vortices and the suitable 
definition of a phase factor whose curl is this approximation of We refer to Section [4] for details and 
simply state the result, again without making the remainder terms precise 



E GP + I TF \ loge| 2 > E GP > E GP + I TF \ loge| 2 + | loge| 



1 



£)TF 2p 



,TF 



(2.29) 



Now an important advantage of our approach becomes apparent, since with such an estimate there is 
no need to rely on compactness arguments to prove vorticity asymptotics. We obtain from (the rigorous 
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version of) ( 2.29 1 an estimate of the norm of the difference /S — /x* and there only remains to use the fact 
that only a small part of the current is removed in the definition (2.24) to estimate the difference between 
fi and jl (see Lemma 5.1| and thus conclude the proof of Theorem 2.2 Note again that, since we expect 
that n contains singularities approaching Dirac masses, the regularization procedure leading to jl seems 
unavoidable in order to properly define the renormalized energy. This justifies the claim in Remark |2|3| 
that ft, — can be estimated in a better norm than /i — fi* . 

To our knowledge this is the first time in the literature that explicit estimates on the rate of convergence 
of a vorticity measure to the minimizer of a renormalized energy are provided. There is the exception of 
[R] but the limit problem was simpler there and its properties not exploited fully as we do here. Common 
to both papers however is the regularization procedure of the vorticity measure that allows to prove 
explicit estimates without using any compactness argument. 
Remark 2.5 (Boundary conditions in (2.6)) 

The reader might wonder why it is natural to use Dirichlct boundary conditions in (2.6), i.e., set h u = 
on dB(R TF ). A first answer is that this is essential to be able to go from (2.25) to (2.26) (see the proof 
of Lemma 5.1). A more physical answer can be provided however: j and j are superfluid currents, so 
they must be thought of as velocity fields. In fact they correspond to phase gradients: where one can 



write u 



ue 



then j 



'Vv? and one should remember that essentially \u\ « 1 except in the small 



region covered by vortex balls. The potential 



boundary condition should rather be thought of as the constraint that 



is therefore defined modulo a constant and the Dirichlet 
hjx is constant on dB(R TF ), i.e., 
on dB(R TF ), where r is the tangent vector. Since j « j should be thought of as l/p TF Vft.^ 
rotated by ir/2 (see (2.25)), this means that we are actually imposing d n ] = on dB(R TF ), a rather 
natural condition if we think of B(R TF ) as the support of the condensate from which the fluid should not 
escape: the superfluid current must be tangent to the boundary of the sample. 

□ 



2.2 Discussion and Extensions 

In many respects, the regime we study here for the GP theory is the analogue of the 'obstacle problem 
regime' of GL theory studied in |SSlj and |SS2[ Chapter 7] (see also |JSlj b Common to both settings 
is the fact that vortices occupy a region whose size grows when increasing O (respectively the external 
magnetic field h ex in GL theory) until it fills the whole sample. The limit problem we obtain is related to 
that of GL theory but has significant differences, mostly due to the inhomogeneous matter density profile 
of GP theory. It is in some sense simpler because it has an explicit solution, but it leads to richer physics: 
the inhomogeneity of the vortex distribution and its progressive homogenization we have discussed before 
are indeed absent in GL theory. 

An analogue of the stability estimate ( 2.28[ ) also holds for the limit problem of GL theory, although it 



does not seem to have been noticed before (see Remark 3[T below). It can be used in the manner we do 



in this paper to obtain slight improvements of the classical results of |SS2[ Chapter 7], like convergence 
of the vorticity measure in better norms and explicit error estimates. 

Our method can also be adapted to treat variations of the physical setting: for instance it is a rather 



simple adaptation to prove the results corresponding to our Theorems 2.1 and 2.2 in the case of an 
annular condensate. Modulo slight modifications one can thus treat the setting of |AAB| in the regime 
where the rotation is above the critical speed for vortex nucleation but of the same order of magnitude. 
Heavier modifications are on the other hand needed in order to extend the results to a regime close 
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to the third critical speecj^jin a flat traps: The GP functional in this case read^] 



S GP m = f dr ( |W| 2 - 2J2**L* + ^L), 
Jb(i) I e J 



(2.30) 



where B(l) is the unit disc, the minimization is performed under a unit mass constraint (that one can 
supplement with a Dirichlet boundary condition) and the regime of interest is Q oc e 2 | loge| 1 with 
e< 1. We refer to |CRY1 [R, CPRYlJ for a more thorough discussion of this model, but only mention 
that if 

fin 2 

e^|loge| Sir 

one reaches a giant vortex phase where the mass is confined to a thin annulus along the boundary dB of 
the trap by centrifugal forces and no vortex is present in the annulus. If 

the present analysis applies and yields the conclusion that vortices are densely packed in an annulus 
included in the bulk of the condensate. This annulus progressively fills the bulk in the limit flo — > 
(that is when fl decreases) while the vortex density, highly inhomogeneous due to a non-constant density 
profile when f^o is not too small, gradually homogenizes. The precise expression of the vortex density 



as a function of the matter density is exactly analogue to (1.15), with the appropriate density p TF and 
cost function F TF (see |CRYj ) . This bridges between [CRY] and |CY[ ICPRYlj where we have proved 
that the vortex distribution is homogeneous when | log e| <C fl <C e _2 |loge| _1 . The borderline case of 
flo = 2(37r)~ 1 (l — o(l)) is considered in [R]. 

As discussed in the introduction of [CRY] , the analysis of the third critical speed in a flat trap is more 
involved than that of the first critical speed treated here. Our method thus needs to be supplemented 
with the tools developed in |CRY|, |R] to adapt to this setting. What makes the application possible is 
that the phase transition happening at the third critical speed can be seen as a fl Cl type transition but 
backwards (vortices disappear when fl is increased). 

Let us also emphasize that the regime corresponding to the third critical speed in a 'soft' trap given 



by a potential such as (1.2 1 is quite different from that in a flat trap (see [CPRY21 ICPRY3j l and shares 
much less features with the fl Cl regime. It is thus unlikely that the methods we develop here can apply 
in this case. 

During the completion of this paper we learned of the recent work [BJOS2] where (among other 
things) the regime we are dealing with has been studied for a 3D condensate. A limit problem is 



derived, formulated in terms of the current (2.14) instead of the vorticity (2.13). The situation seems 



more complicated in 3D, and nothing as explicit as formula (1.15) appears to be derivable from this 
limit problem. The fact that, when fl 3> |loge|, the vortex distribution becomes homogeneous (in fact 
constituted of many densely packed and uniformly distributed straight vortex lines parallel to the axis of 
rotation) is proved however, confirming our results in the regime where the rotation speed largely exceeds 
the critical one. Whether our approach, supplemented with the tools of [BJOS1 , can be generalized to 
three space dimensions and complete the results of |BJOS2j remains a question for future investigations. 

We finally remark that vortex patterns inhomogeneities play a crucial role in a very different regime 
than that under consideration here, namely in the vicinity of the maximum rotation speed attainable 
in a condensate confined by a purely quadratic trap |AB1 IABD] . In that case the inhomogeneity only 
manifests itself close to the boundary of the condensate and is thus hardly observable, but it has the 
important effect of modifying the density profile of the fluid. 



6 For a description of the physics of the three critical speeds of GP theory, the reader may want to refer to |CPRYl1 
ICPRY2I . 

7 Note the slightly different units as in CRY, CPRYl], i.e., mass equal to 1/2 and angular velocity 2f2. 
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2.3 Plan of the Paper 

The rest of the paper presents the proofs of our main results. We start by analyzing in Section [3] the 
properties of the renormalized energy. Section [4] contains the construction of our trial state and the 
evaluation of its energy. The general technique is not new but, since our method allows to deduce from 
energy bounds a quantitative estimate of the rate of convergence of p, to we make an effort to obtain 
precise estimates of the remainders in the energy upper bound. Section [5] is then concerned with the 
energy lower bound and the proof of Theorem 2.2 An appendix gathers technical results used in several 
places of the proofs. 



3 The Renormalized Energy 



In this section we focus on the study of the renormalized energy of vortices. For further convenience we 
define the limit functional for a slightly larger class of densities p and measures v as 



V 



1 

2f> 



1 



V/i„| 2 + -p\ v \+Fv 



where 



K = 0, 



in V, 
on dT>. 



(3.1) 



(3.2) 



We will not strive for the most general assumptions allowing the study of such an energy functional, but 
instead state a theorem that is sufficient for the purpose of proving our main results. From a mathematical 
point of view, this means that the functions p and F and the domain T> appearing in the definition above 
are not necessarily the ones that appear in the preceding sections, but the assumptions we impose on 
them are inherited from the physical features of our original GP theory. Note however that we do not 
assume in this section that the problem is radial. We do not stress the dependence of I p on F because 



in the applications we have in mind F is related to p as in (1.10). 



The following theorem states the existence and uniqueness of the minimizer of the functional I p in 
the natural energy space. These are somewhat classical results reminiscent of potential theory |ST|IStaj . 
more important to us are the explicit formula (3.6) for the minimizer and the stability property (3.8) 



that is the key input in the proof of Theorem |2.2| 

Theorem 3.1 (Minimization of the renormalized energy). 

Let V be a regular open subset ofM 2 , p E C 2 (2?), p > 0, F g C 2 (T>). We assume that 



1 



-p-F>0, inV, 



(3.3) 



and define 



H p :=-p + F 



(3.4) 



1. Existence and uniqueness of a minimizer. The functional X p has a unique minimizer u p in 
the class of measures 



J( p (p) 
It is given by the formula 



1 



{H„<0}, 



(3.5) 



(3.6) 
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and the ground state energy is 



I p = l p [p p ] = - ( Hpfip. 
2. Stability of the minimizer. For any v 6 j$ p (T>) 



(3.7) 



v 



2p 



\Vh 



(3.8) 



The main physically relevant assumption we make here is (3.3 1, which indicates that negative degree 
vortices are not energetically favored. In our setting it will always be satisfied because typically p > 
and F < with equality only at the boundary of the domain. One may certainly prove a related theorem 
in the case where (3.3) does not hold but it is not our concern here. 

Proof. We split the proof in four steps. 

Step 1 (Existence and Euler- Lagrange equation). It is not difficult to prove the existence of a minimizer, 
we will thus skip most of this discussion. Let us just note that with our assumptions 



2pM> / ~|v/g 2 + c'iHff- 1 

JV 4 P 



(25) 



C 2 IMI(C c i(X>))* + / P\ v l 



V 



with C x = (2 

sup r gp p) 1 and C2 — sup r g.p (|p — F). Recalling the embedding of (C* (£>))* in H~ l {V) 
it is then easy to deduce bounds on the minimizing sequences and conclude by lower semi-continuity 
arguments. 

Considering now a variation of the form (1 + tf)p p , f £ C (T>), and noticing that (3.2) implies 

1, 



v P 



huV, 



we see that the Euler-Lagrange equation of the minimization problem takes the form 



for any / e C°{V). Writing 

we deduce from the above that 
-H, 



h P P P P + \p\Pp\ + F Pp ) / = 



p p = fj+ - p p , with p+, p p > 



(3.9) 



(3.10) 



p p ---p, on supp^/ij J, 
Mp = |/ , -- F ! on supp (^"), 



Pp = V y^Hp) > on supp (pj 



Pp 



P V(\P-F) 



(3.11) 



on supp I n 



Note that this implies in particular supp (p^) C |v ViJ p ) > 0}. 

Step 2 (Explicit formula for the minimizer). We first prove that p~ = by using the trial state 



Atrial — Pp 



Uniqueness of the solution to (3.2) implies that 



Hp, on supp (pj) 
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Atrial being simply extended to the whole domain D by requiring that V(p 1 V/i Ptrial ) = on T> \ 
supp (/il"). Recalling (3.9) and (3.11), we thus have 



p [/^trialj 



whereas 



Using assumption (3.3), this clearly means that it must be p p = 0. The minimizer p p > is then 
completely determined by the set supp (p p ) C {V(^VHp) > 0} and we have (recall (3.11)) 

1 Jsupp(p p ) VP / 

from which it is easy to deduce that, in order for \x p to minimize I p , it must be 

supp (n P ) = (v (-VH P ) > o 1 n {H p < 0} , 



and (3.6) follows. 

Step 3 (Consequences) . We first claim that 



hn P > ~H P , a.e. in V. 



(3.12) 



Indeed, we have equality in supp(/i p ) = {V(^Vi?p) > 0} n {H p < 0}. Let us now suppose that ( |3.12 1 
does not hold in some sub-region of supp(/i p ) c and pick some v > with support in this region. Consider 
the trial state 

Atrial = Hp + tv 

for t > 0. Since p p and v have disjoint supports 

l/Viall = lMpl + 

Also we remark that for any pair of measures p.\, P2 € ^{T>) 

/ -V/i Ml • V/i P2 = / MiV*' 
JX> P Jv 

which allows to compute the energy of ptrmi arid find 

I P Atrial] =X P K]+i 2 / ^K\ 2 + tf (H, 
Jv Z P Jv 



(3.13) 



Since the last term is negative by assumption, for t small enough we would obtain I p [/itriai] < 1p[pp\, 
which contradicts the fact that p p minimizes I p . 
Next we prove that 

-p-F-h Pp >0, a.e. in supp(^ p ). (3.14) 

The argument is again by contradiction. Suppose that in some subregion of supp(/i p ), we have p — 2F — 
2h Pp < and pick some positive measure v € ^(D) whose support is included in this region. Consider 
the trial state 

Atrial = Pp — tV 
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for some t > 0. We have 



\fi p - tu\ < \p p \ + tv. 



because of (3.2). Computing Xp[/itrial] by expanding the quadratic term and using (3.131, we thus find 



Atrial] <l p [fl p ]+t 2 J^^-\VK\ 2 +t Q 



p-F-K 



from which we immediately deduce that, for small enough t, we would obtain X p [//trial] < ^pIMpL which 
is a contradiction. We conclude that (3.14) must hold true. 

Step 4 (Stability) . We are now ready to prove the most important point of our theorem, namely the 
stability property (3.8). We write any v g jM{T)) as 

V = Hp + Vi + l>2, 

where supp(^i) C supp(/i p ) and supp^) C supp(p p ) c and note that 

\u\ > Hp + v\ + \v 2 \. 



Since p p > 0, using (3.13) again, we get 

!>] > 1 P W] + j- p \Vh Ul+ „ 2 1 2 + j Qp + F + hA v x + J \p\v 2 \ +f v (F + K P ) "2, (3-15) 



but p + 2F + 2h P] = on supp(^i) by (3.11 ). Thus, decoupling v-i into its positive and negative parts: 



za> = v 2 



we have 



lp[v}>lp[p p ] + jT i|VA 



v-r+v-z I 



with 1*2 , v 2 > 



-p + F + hp p 



-p-F-h 



The last two terms of the above expression are positive because of (3.12) and (3.14) and this leads to 

Q . □ 

Remark 3.1 (The obstacle problem of Ginzburg-Landau theory) 

As already mentioned, the problem studied in this section has a clear connection with the limit problem 
obtained in the study of type-II superconductors in the first critical field regime [SSI, SS2 , which in 
turn is connected to obstacle problems |BSj . The energy of a vortex density in a type-II superconductor 
occupying the 2D domain T> can be approximated in suitable units and variables by 



rGL 



where h Ll solves 



and A > is a parameter. 



-Ahp + hp 



p, in T>, 
on &D, 



(3.16) 



(3.17) 



It is known |BS| ISSl] that (3.16) has a unique minimizer p GL , expressed in terms of the solution of 
an obstacle problem, which is a particular type of free-boundary problem. As far as we know, no explicit 
formula for p GL exists for generic domains T>: it is constant in a subdomain of the sample T> and zero in 
the rest of the domain but the boundary between the two regions is not known explicitly, although much 
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can be proved about it (see }SS2l Chapter 7] and references therein). In the special case of a spherically 
symmetric domain, one can show by uniqueness of the minimizer that p Gh is in fact radial and therefore 
almost explicit. 

Using a method similar to that we used for the proof of Theor em|3.1[ one can show the following: for 
any v such that /i GL + v is in the energy space corresponding to (3.161, 



1 GL [p GL + v] > Z GL [M GL ] + \ J (WK 



\h u -l\ 



(3.18) 



which is the equivalent of the stability estimate (3.8). A convenient way of seeing that (3.18) holds true 
is the change of variables p — > p, — 1, — > — 1, which yields as equivalent problem the minimization 
of the functional 



l[p} = 



1 



v 



|m + i| 



\Vh„ 



l 



V 



|m + i| 



+ h^p 



where now 



-Ahn + hu = p, in T>, 



on dV. 



h M = 0, 



This form has the advantage of being closer to (3.1) and facilitating the proof of (3.18), but the physical 



interpretation of the problem is more transparent in (3.16) 



□ 

Now we apply the result proven above to the functional I TF defined in (2.4). In this case both F and 
Hp can be explicitly computed: 



F (r) 



dttp r *(t) 



R 



TF' 



(i? Tp2 - r 2 ) 



i? TFS+2 - r s + 2 



(3.19) 



H"(r) = -p"(r) + F"(r) 



R 



■'J L' ,s ^TF 2 _ r 2^j _ 1 ^ 



R 



TF 



s+2 



n s+2 



(3.20) 



There is a certain freedom in the choice of the integration domain V (see Remark 3J2 ) in X TF [v] but for 
clarity we set V = B(R TF ). 

Corollary 3.1 (Minimization of I TF [i'}). 

Let I TF [i'] be the renormalized energy defined in (2.4) with p TF given by (1.6). There exists a unique 



minimizer p± in the class of measures (3.5), which is radial and absolutely continuous w.r.t. the Lebesgue 

(3.21) 



measure, i.e., there exists a continuous radial function m*(r) such that 

p* = [m*(r)] + l {ffT F< 0} dr 
and 



m*(r) = ^ r 2 log(p TF (r))+2Q 
If in addition Hq > the support of p* satisfies 

? supp(^) = S(i^) C 6(i? TF ) 

for some < i?* < i? TF swc/i that 



R±R 



TF 



(3.22) 

(3.23) 
(3.24) 



when S7q ~^ 00 • 
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Remark 3.2 {Domain o/I TF [i/]) 

A simple inspection of the proof of the above results shows that the minimization and thus the ground 
state energy and minimizer are to some extent independent of the domain T> in the definition of I TF \v\ . 
More precisely one obtains the same minimizing /i* and the same minimum energy for any T> such that 

B(R+) cV c B(R TF ). (3.25) 

Indeed /i* is a minimizer as long as supp(/z*) C V, so that ( |3.23 1 yields the condition (3.25). 

□ 



Proof. In order to apply Theorem 3.1 we can take any T> C B{R ), since by construction F(r) < for 
any r < R TF and F TF (R TF ) = 0. Moreover p TF and F TF clearly satisfy all the requirements of Theorem 
|3 - 1 1 and the first part of the statement is thus proven. 



The explicit expression of /i* is provided by (3.6) and since both p TF and H TF are radial, /x* must be 
so as well. A simple computation yields to*, which is a well-defined continuous function for any r < R TF . 
The condition _ff TF < guarantees that the support of fi± is strictly contained in B(R TF ): Indeed one 
has# TF (i? TF )=0but 

H TF '(R TF ) = lp TF '(R TF ) = —sR™ 8 - 1 < 0, 

which implies that there exists some Ri < R TF , such that H TF (Ri) = and H TF (r) > for any 
Ri < r < R TF . The fact that R± > follows from the analysis of H TF at the origin: using (1.7) and 
(1.13), we compute 

H TF (0) = -j A TF - -^fi = (O x - fi ) < 0, (3.26) 

4 Z7T ZTT 

by assumption. Moreover a simple analysis of the derivative of H TF , i.e., 

shows that H TF (r) = for r = 0, R for some R > 0. It is positive between and R and negative for 
R < r < R TF . In particular this, together with the discussion above, implies that H TF has a unique 
maximum in [0, R TF ] and, since H TF is negative at the origin and vanishes at the outer boundary, one 
can conclude that 

H TF (i?!)=0, H TF (r)<0, for any < r < R v (3.27) 
To complete the proof it suffices therefore to exploit the explicit expression of to*, i.e., 



r s-2 ( r s + s(s _ 1)i? TF^ 

8p TF (r) 2 



m*(r) = V T w_, a L + 2n , (3.28) 



which easily implies that to* is decreasing since the first term on the r.h.s. of the above expression is. 
Moreover 

to*(0) = fl > 0, m*(r) > -oo, (3.29) 

so that the subregion where to* is positive is another ball of radius i?2 < i? TF and 

R+ = mw[R 1 ,R 2 ]. (3.30) 

To see that i?* — > i? TF in the limit f2 — > oo, it is sufficient to note that the maximum point R of H TF 
converges to i? TF as fto — > oo, which implies that sup 0<r<fl TF H TF — > in the same limit. Hence 

Ri > R TF (3.31) 
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and an analogous statement holds true for R 2 , i.e., 



R 2 > R 

Sin— ^oo 



TF 



(3.32) 



as it immediately follows by noticing the first terms in (3.221 is independent of Sl , while the second is 
linear in Qq. □ 

A simple computation yields the equations satisfied by R\ and R2: setting r.; :— Ri/R TF , i = 1, 2 (so 
that < ri < 1), one has 



l-n 



1 - r 2/s - 

1 ' 1 s+2 



n R 



1 (»-2)/» s(s-l) + r 2 

2 2 (l-r a )2 



n R TF ' 



(3.33) 



which allows to obtain a relation between r\ and r 2 independent of fio- The computation of T\ or r 2 or 
their ratio can not be made explicitly but the above equations could be tested numerically to study the 
dependence of the support of /i* on the parameter s, that is whether R+ is given by R\ or R 2 . In the 



harmonic case s = 2 however both equations (3.331 are solvable and one obtains 



Ri 



n R TF ' 



<2 



R2 



R TF \ An Q R^ 2 




2^ i? TF " 



(3.34) 



Note that, by taking the limit il — > 00 of the above expressions, one easily recovers (3.24 1. 



4 Energy Upper Bound 

The main result of this section is stated in the following Proposition |4. 1| and is the proof of an appropriate 
upper bound for the GP ground state energy. 

The result is obtained by testing the GP functional on some explicit trial function Atrial- The most 
difficult part is the construction of such a function because of several requirements it has to fulfill in order 
that its energy is suitably close to the ground state energy. There are indeed two main constraints: On 
the one hand the modulus of Atrial must be approximately equal to g, in order to recover the leading order 
contribution E GP , but, on the other hand, its phase has to contain a large number (of order 0( \ logeQ) 
of vortices of unit degree, distributed according to the minimizing measure /i* given by (3.6) on a scale 
I log e I . In addition Atrial must belong to ^ GP , which in particular implies that it must be normalized in 
L 2 (R 2 ). The rest of the proof is just the evaluation of £ GP on Atrial- 

Proposition 4.1 (GP energy upper bound). 

If VL = fiol logej, with Hq > fl± as e — > 0, then 



< E 



GP 



I T¥ \\oge\* (1 + 



log I log £ I 

I logerl 1 / 2 



(4.1) 



4.1 The Trial Function 

As we have anticipated at the beginning of the section, the most difficult part of the proof is the identifi- 
cation of the proper trial function to use in the estimate of the ground state energy. In order to simplify 
the analysis however we first extract from E GF the leading order term, i.e., the energy E GP . Like in the 
lower bound, this can be done by a splitting trick, i.e., setting , 3> GP =: gv, where g is the vortex-free 
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profile given by the minimization (2.3) and v some locally bounded function. Exploiting the variational 



equation for g, one obtains the identity 



E GP = E GP + £[v], £[v] := J dr { \g 2 \Vv\ 2 - % 2 r^ • (iv, Vv) + e~ V (l - |v| 2 ) 2 } , (4.2) 



where we recall that 

i , , ,. _ , i 



r := (—y, x), (iv, Vu) :— - (vVv* — v* Vu) . 

It is clear that in order to obtain an upper bound for E GP , it is sufficient to provide a suitable trial 
function ftrial to evaluate £[v]. Note that v inherits an L 2 -normalization from ^> GP : 



/ dr g 2 \v 



2|„.|2 _ i 

r , 



and the same must apply to "trial- We thus set 

"trial (r) := c £ £(r) exp {Atrial} , (4.3) 

where £ is a cut-off function, Atrial is a phase factor we are going to define and c e is a normalization 
constant. 

The phase factor of "trial must contain a number O(|loge|) of vortices of unit degree distributed 
according to the measure loge|. For further convenience we start by defining the phase inside a ball 
B(R < ) of radius 

R K :=i? TF -e 2 / 3 |log £ | 2 / 3 , (4.4) 
which is strictly contained inside the TF support B(R TF ) and has the crucial property that inside B(R 



the pointwise estimate (A. 5 1 applies and g 2 ~ p TF : There we define Atrial as the solution of the equation 

N c 

V^rial = -P TF V x /i^ rlal) forreB(R < )\\jB( ai ,e), (4.5) 

i=l 

where the potential h Vuiai solves the differential equation (see below for further details on the existence 
and properties of such a solution) 

-V (p TF_1 V^ trlal ) = Atrial, in B(R < ), 

-Ah UtM =0, in B(R > )\B(R < ), (4-6) 

^triai = 0) 011 9B(R > ), 

for some R > > R TF given below. Note that we modify the equation to A/i„ trial = in the boundary layer 
B(R>) \ B(R < ). This is to avoid problems due to the vanishing of p TF on dB(R TF ). Other strategies are 
possible, see, e.g., [B] Section 2]. The outer radius i?> is larger but close enough to R TF : for technical 
reasons which will be clearer later and are mostly related to the estimate of the interaction term and the 



exponential decay (A. 8), we pick 

R> :=i? TF + e 2 / 3 |log £ | 4 / 3 . (4.7) 

The function Atrial is a smooth approximation of a measure given by the sum of Dirac masses at the 
points a 4 , i = l,...,N e , e.g., 

2 Nc 

Atrial := ~2 le ( a - £ )' ( 4 - 8 ) 

i=l 
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i.e., SLi are the positions of vortices and N £ their total number. We have also implicitly assumed that 
inf |aj — a.j\ > 2s and < R < — e (see below for further details on the distribution of points). The crucial 
property of Atrial is that, given any contour C C B(R < ) encircling only one ball B(aj,e), one has 



do- d T 



'trial 



da d T 



'trial 



3B(a,,e) 



JdB(a.i,e) 



8 h 



[ drv(p TF Vfc^J = / Vtxial 



2tt, (4.9) 



B(a i)£ ) 



since the integral of curl(V^triai) over the area delimited by C and dB{&i,e) vanishes thanks to (4.6). 
Hence Atrial is a well-defined phase factor for any r G B{R < ) \ U£>(a,;, e), whereas inside the vortex balls 
Atrial is not the phase of any function and we will have to use a cut-off function £ to exclude that region. 
We need however to define the trial function and therefore its phase up to R > because otherwise the last 
term in £ [v] (see (4.2)) would give a too large contribution. For the phase in the boundary layer we set 

V^triai = -p TF_1 (i?<)V ± /v trial , for r G B(iZ>) \ B(R K ). (4.10) 

In fact since the boundary layer B(R>) \B(R < ) is vortex- free by construction, it is not difficult to realize 
that the differential equation (4.6) can be solved explicitly there and 

friw = clog(r/i?>) , 

where the coefficient c is fixed by imposing continuity of d n h Ui , tial on dB(R < ). On the other hand if wc 
integrate (4.6) in B(R < ) we obtain 



/ dr V (p TF 1 Vh v 



TF" 



(R< 



dcr d n h Vtiial 



dB(R < 



which combined with the explicit expression of h UtiiBl for r > i?< yields 

K tI1 Jr) = - N e p TF (R^ log (r/Ry), for r G B(R>) \ B(R<), 

and therefore 



V0 



'trial 



= N F _ 



for r G B{R > )\B{R < ). 



(4.11) 



(4.12) 



(4.13) 



One can verify that this extension defines a phase in B(R > ) \ UB(&i,e): the property (4.9) is still valid 
and, for any contour C C B(R>) containing all the vortex balls B(a.i,e), one has 



dcr d T Atrial 



dcr d T c 



'trial 



(4.14) 



ae(i?<) 



thanks to (4.13) and the fact that the r.h.s. of (4.5) or (4.10) is irrotational outside C. 

We can now turn to the cut-off function £. Its role is twofold: on the one hand it is needed to restrict 
the integration domain in £[v] to S(i?>) and, on the other, regularizes Atrial inside the vortex balls to 
ensure that Vtriai is a one-valued function. Concretely we set 



£(r) :=fo(r)n£v(|r-ai|), 



(4.15) 



;=i 



with £o and £ v smooth functions satisfying 



&(r) - 



for r > i?>, 
for r < R y — : 



2 / 3 I 10' 



0, for r < e, 

1, for r > 2e, 



(4.16) 
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and such that their gradients satisfy the conditions 

HVfolL <0(e- 2/3 |loger/ 3 ), 



|Ve v |L < Oie- 1 ). 



(4.17) 



With these definitions we can estimate the normalization constant without saying anything more 
about the distribution of points but only assuming that N e — C(|loge|): since the cut-off function 
satisfies £ < 1, one clearly has c e > 1 and 



dr g 2 tl 2 > c 2 / dr g 2 

1 JB(fl>-e 2 / 3 |loge|- 1 )\uB(a l ,2e) 



>C 



1 



'\B(H>- E V8|lo ge |-i) 



dr 3 2 -C(e 2 |loge| 



>c 2 [l-0( £ 2 |log £ | 



thanks to boundedness and monotonicity of g and the exponential decay (A.8), yielding 

,2/ 



g 2 (r) = Oie 00 ), uniformly in r > R TF + 0{e 2/3 \ loge| 4/3 ). 



Therefore we conclude 



I<c 2 <l + O( e 2 |loge|). 



(4.18) 



(4.19) 



Finally we discuss the distribution of vortex points : the final goal is the reproduction of the density 
provided by the minimizing measure loge|. We proceed as follows: for any k € N larger than some 
given ko > 0, we set 

Qk ■= -J=t, (4-20) 

vl lo s £ l 

and denote by Ck and Nh the circle dB(gk) and the number of equidistributcd points we will put on Ck 
respectively. Recalling (3.22), we put uniformly distributed points on Ck in such a way that the total 



number is (with |_ • J standing for the integer part) 



N, 



2tiV| loge|£> fe m*(£> fe ) 



(4.21) 



provided Ck C supp(^^).We denote by K. C N the set of integers k such that Ck C supp(/i*), so that one 

clearly has 

sup k < R^\\ogE\, sup |a Jfc | < R+, (4.22) 

keK. k£K 



thanks to (3.23). The lower bound k > k is needed in order to ensure that on the smallest circle there 



are sufficiently many points (e.g., Nk > 4 for any k > ko). The point distribution naturally provides a 
decomposition of Ck into sectors di k < d < t?i fc +i, ifc = 1, • • • , Nk, centered at the points a.i k and of equal 
arc length given by \Ck\/Nk, i.e., for k e JC, 



0; 



2ir/Nk 



O fc := ih k+ i 

The total number of points is of order 0(\ logsQ: 

7V £ = ^iV fe = 0(|loge|). 



(4.23) 



(4.24) 



keK 
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Moreover z^riai satisfies the required property, namely i-'triai — |loge|/i*. More precisely for any radial 
test function $(r) € C 1 (S(i? TF )), 



/ dr $(r)v tlial = (1 + 0(e)) ]T ]T $ (a, J = (1 + 0(e)) £ iV fc $ ( fffc ) 

= 2tt (1 + 0(e)) VUoge| ^ e k m*(e*)* (g*) + (V|loge|) 
= \loge\ [ dr$(r)m*(r) + 0(V|loge|) = |loge| / + (Vflogt 

■/suppfuj \ 



where we have used the one-dimensional Riemann sum approximation, 



2^| log er 1 / 2 ^^/^)- f dr 



fir) 



<C||/ / (r)|| 00 |lo ge |- 1 /2 j 



the fact that by hypothesis ||V$|| < C and the gradient estimate 



1 



1 



ll m i(0llL-( B upp(/O) " 9 SUP l^rlogp | = " SUp 



TF2 TF'» , / TF'\ 3 o TF TF' TF" 



by (3.23) 



4.2 Evaluation of the reduced energy £[vtriai] 

Now we are going to estimate each term of the energy E [vtriai] separately. 
We start by considering the kinetic term: Noticing that one has 

|V<w| 2 = c 2 |V£| 2 + c 2 e e |V0 trial | 2 < (1 + 0(e 2 | loge|)) (|V£| 2 + £ 2 |V0 tria i| 2 ) 



by (4.19), we obtain the bound 



/ dr g 2 |V Wtria i| 2 < (l + 0(e 2 |loge|)) [ / dr.g 2 |Ve| 2 +/ dr«? 2 |V0 tr 

JM 2 Ur 2 JB(R > )\uB(a,,s) 



Lai 



(4.25) 



< C, (4.26) 



< (l + 0(e 2 |loge|)) / dr 9 V F 2 |V/i„ tr ,J 2 + 0iV 2 / dr r~ 2 g 2 + 0(| loge| 

JS(fl<) iB(H > )\B(fl < ) 



< (l + 0(e 2 / 3 |loge| 2 / 3 )) / drp TF 1 | W rial | 2 + 0(| loge|), (4.27) 



where we have used the upper bounds (4.17), the estimate of the total number of points (4.24) and (A. 5 1 



for the term located in B{R < ). In the boundary layer R < < r < i?> we have exploited the explicit 
expression for V Atria l (|4.13 l and the pointwise estimate (A.8). We estimate now the interaction term, 
i.e. the last term in (14. 21): 



e- 2 / dr/fl-h 



trial | 



< £ 



! \B(i?>- £ 2 / 3 |loge|-i) 



dr g A + £- 2 



B(i? 5 



dvg ( 1 — |v t rial| 



< C*e~ 2 (c 2 - if + 0(e°°) < 0(e 2 | loge| 2 ), (4.28) 
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thanks to the exponential decay of g ( |A.8 ) and (4.191. 

The angular momentum term (second term in (4.2 )) can be integrated by parts: 



Q / dr g 2 r^ ■ (iv tlia i, Vu t riai) = I loge| / dr F g (r) curl (w tria i, V 
Jr 2 Jbir-,) 



where 



F g (r) := -Q / dttg 2 (t) 



Wtrial) , (4.29) 



(4.30) 



Note that in the integration by parts the boundary term vanishes thanks to the fact that F g vanishes 



there. In fact in the boundary layer £?(i?>) \ B(R < ) it is more convenient to integrate the r.h.s. of (4.29 1 
by parts back: 



loge| / dr F s (r)curl(w trial , V^triai) 

Je(i?>)\8(fl<) 



-fill log e| 



B(it>)\e(« 



< £ 2/3 



e(i?>)\e(fl 



dr g 2 Y ± ■ (Wtrial, VDtrial) ~ | log e\ I d(T F g {R < ) (Atrial, &r«tiial) 

JdB(R^) 

dr 3 2 |V0 tI , al | 2 + Ce- 2 / 3 |loge| 2 / 

<) JB(fl>) 

< £ 2 / 3 / dr.g 2 |V0 tl , al | 2 + O( £ 2 / 3 |lo ge | 10 / 3 ) (4.31) 



drr 2 ff 2 + |F s ( J R < )|0(|log£|) 



by Cauchy-Schwarz inequality and the pointwise estimate (A.8|. The same estimate also implies that 
\F g (R < )\ < Ce 4 / 3 | loge| 4 / 3 , which has been used above as well as the fact that 



d<J (Atrial, 9 T Utrial) = C\ 



dB(R < 



[ da 5 T tl ,ai = 2^c 2 deg {0 tl ,ai, 5B(i?<)} = 0(| log e|), (4.32) 
Jae(fl<) 



by (4.141 and ( 4.24 1 . The r.h.s. of (4.31) can then be safely incorporated in the remainder in (4.27) and 



it is sufficient to consider the r.h.s. of (4.29) only inside B(R < ), where £o = 1 an d 



curl (retrial, Vvtriai) = c 2 curl (V(/> trial 

) -0, in |JS( ai ,2e) f)6(i?<), 



(4.33) 



so that wc have 



loge| / dr F s (r)curl (TO tri ai, V^txiai) = | log e| / dr F g (r) curl(wt ria i, Vu t riai) 

Je(fl<) ^UB(aj,2e) 

= (1 + 0(e)) | log e| Y»i) / dr curl (£ 2 V0 tria i) 

i=1 Je(a,,2e) 

= (1 + 0(e)) | log £ | y2 F M / da Atrial < 2tt| loge| V ^(a,) + 0(e| loge| 2 ), (4.34) 



i=l 



where we have used the bound ||VF ff ||oo < C to estimate 



sup |J7 9 (r)-F ff (a0|<Ce|.F fl (ai)|. 

reS(a;,2e) 



(4.35) 
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since 



F g (a t )\>cJ^ < &ttg 2 {t) > C (l - C(e 2/3 | log e| 2 / 3 )) J < dt p TF (t) > C > 0, (4.36) 



thanks to (3.23), (4.22) and (A. 5). The replacement of F g by F can be done in a similar way: 
iF^-F^ia^^C f R> dtt \g 2 (t)-p TF (t)\ < 0( £ 2 / 3 | log £ | 2 / 3 ), 



(4.37) 



by (4.22) and ( A.5[ ), so that using (4.36) one obtains 

\F g (aA-F TF (a t )\ <0(e 2 ^\\oge\^ 3 )F TF ( ai ), for any i = 1, . . . , N e 
Combining ( |4.27| and ( |4.28| ) with ( |4.29[ ), ( |4.31[ ) and ( |4.34[ ), we obtain 

1 



£ [atrial] < l + 0(£ 2/3 |loge| 2 / 3 ) 



V^ triaI | 2 + 2n\ loge| £ F 3 (a,) + 0(| loge|), (4.38) 



O TF 'Atrial 

0(ii<) 2 P 

so that in order to complete the proof it suffices to estimate the first term on r.h.s. of the expression 
above. To this purpose we first note that the potential 



K tM := K %M + N eP lb (i?<) log (i?</i?>) 



solves 



( TF — 1 
[P 



Atrial = > 



V tri al, in^(i?<), 

on dB{R K ), 



(4.39) 
(4.40) 



since h Vtlial is continuous in the whole of B{R > ) and for r > i?< it has the explicit expression (4.f2). 
Obviously one also has 



dr -4p |V^ trial | 2 



dr -=|V/i 



B(fl<) P 



XF I "trial I 



It is now convenient to introduce the Green function 

ty(x, y ):=A- 1 (x, y ), A:=-v[p TF ^ 



(4.41) 



i.e., the integral kernel of the inverse of the differential operator A. More precisely Q solves the differential 
equation (we postpone to Lemma 4.1 below the proof of the required properties of Q) 



P 1 * (x)V x £(x,y) =<J(x-y), for x, y e B(i?<), 



^(x,y) = 0, 



for x <= dB{R<). 



(4.42) 



By means of G, we can rewrite the first term on the r.h.s. of (4.38) as 
1 



B(R<) P 



TF l^^triall 



0(x, y)ftrial(x)^trial(y) 



B(fl<) ./B(fl<) 



^/ dx/ dy^x,y) + 4E/ dx / «Jy5(x,y). (4.43) 

£ i=1 -'B(a l ,£) JB(&i,E) £ JB(tu.e) JB(sn.e) 
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For the estimate of the diagonal term we can use the inequality (4.53) proven in Lemma 4.1 below, 
which yields 



E/ dx / dyg(x,y) <-^X] / dx f dyp TF (y)log|x- y |+0(|log £ |) 
<-(l + 0( e ))_ ^p TF (aO / dx / dy log|x-y|+0(|loge|) 

< 2^]T /9 TF (a l )|log£| +0(|loge|), (4.44) 



where we have used the properties of harmonic functions and in particular the mean value theorem to 
compute 



/ dx / dy log |x - y| 

1 



2 7o 



AO 



dd' \ dgg / dg' g' log g 2 + g' - 2gg' cos (tf - ti>) 



= tt e \ log el 



4ir 2 / dgg / dpY log Q + / d^V log g' 

J0 Uo Jg 



1 



7r 2 e 4 . (4.45) 



To estimate of the difference p TF (r) — p TF (a,i), we have used the bounds || Vp TF || 00 < C and 

p TF {a t ) > C > 0, for any i = 1, ... , 7V e , 



(4.46) 



which is a simple consequence of the explicit expression of p , (3.23) and (4.22). 



The off-diagonal term on the other hand can be dealt with by using the inequality (4.54) proven 



in Lemma |4.1| and a Riemann sum approximation. The distribution of points {a^} identifies a cell 
decomposition of S(i?*): one can associate with each sector a two-dimensional cell Qi k as 



with volume 



(r,0) : J? G 



,r G 



(A-^llogel-Va.^ + ijIlogel-Va]!, 



ICiJ = 



> 



x/Uoge] AffeVTlogE] mJe fe )|loge|' 



(4.47) 
(4.48) 



by (4.20), (4.21) and (4.23). Given a it will be convenient to denote n(ik) the set of indices jh such 



that Qj h is a nearest neighbor of Qi k in the cell decomposition. Indeed, <?(x, y) is singular for small 
|x — y | and it is therefore useful to distinguish those pairs of cells Qi k , Qj h for which the distance |x — y | 
can be arbitrarily small for x G Q«fe,y G Qj h . 



We first write, using (4.54) and the fact that the distance between and a.j is at least of order 
loge| -1 / 2 for i ^ j, 



V/ dx / dyS(x,y) = (2^) 2 V S(a l5 a,) + 0(e| log £ | 2 log | loge 



e i~xJB( ai , s ) 



(4.49) 
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Then we neglect the neighboring points in the sum, using (4.53) and the fact that ja^ — a Jh | is of order 
| loge|~ 1/2 for j h e n{ik): 



N c 



N k N h 



J2G(^^)= E E E 0(m t) a,J + 0(|loge|log|loge| 



'■j=i 



k,h£Ki k = t jh = i 



Next we use (4.541 to write (here its is crucial to have excluded from the sum the nearest neighbors of 
each point to have lower bounds on |x — y|) 



( 2 -) 2 E E E £K,a,j= E wrk-\Y. E / dx / d ^y) 

k,heKi k = l j h =l k,h£K 1 k 1 1 Jh 1 i k = l j h = l J Qi k J &ih 



2 N k N h 



k,h£Ki k = l ifc=l 



O(|loge| 3/2 log|loge|). (4.50) 



Now, using (4.53) again, one can see that neighboring cells do not weigh too much in the integrals: 
E fTTTTTj I E E / dx / dyS(x,y) = 0(|lo g£ |log|log E |), 

which allows to conclude, gathering the above equations and using a Riemann sum approximation again, 

N e 



E/ dx / d y g(x, y ) = |loge| 2 f f £(x, 



0(|log £ | 3 / 2 log|log £ |). (4.51) 



We can now go back to (4.38) and replace (4.44) and (4.49), finally obtaining 

N e r. 



£ Atrial] < 27r|loge|^ 



i=l 



1 



JV e 



+ (2tt) 2 £ Q(ai,aj)+0(\loge\) 



< 



2^|loge| Vi/ TF ( ai ) + |log £ | 2 / / g(x,y)^(x)^( 2 /) + 0(|log £ | 3 / 2 log|log£|) 
<|log £ | 2 / dr (4^|VV| 2 + i? TF (r)m4r)l+0(|log £ | 3 / 2 log|log £ |) 

JB(fl TF ) IP J 



/ TF | log £ | 2 + C(| log £ | 3/2 log I loge|), (4.52) 



4.1 



is 



where we have used (4.25) for the term involving H TF . Recalling (4.2), the proof of Proposition 
now complete. 

We conclude this section by proving some technical results about the Green function 0(x, y) used in 
the proof above. We recall that Q is defined as the solution of (see (4.42)) 



-V, 



„TF" 



L (x)V x g(x,y)j = *(x-y), for x, y e B(R<), 
i 0(x,y)=O, ' forxe dB(R<). 
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Lemma 4.1 (Properties of the Green function Q(x,y)). 

There exists a unique solution C/(x,y) to (4.42) in W 1,P {B{R < )®B{R < )) for any 1 < p < 2. It is positive 
and symmetric and 

1 

2tt' 



sup 

x,yeB(-R< 



5(x,y) + ; fp TF ( 2/ )log|x-y| 



< c 



(4.53) 



for some fixed constant C > 0. Moreover, for any i,j € N and y 6 Q, 



xGQ, 



sup |V a £(x,y)| < C sup |log|x-y|| +1) , 



x6Q; 



(4.54) 



where the cell Qi is defined in (4.47|. 

Proof. Most of the properties of Q follow from standard arguments applied to the analysis of the elliptic 
operator A (see, e.g., [Staj ). Moreover analogous results including (4.53) have already been proven in [R] 
Lemma 2.4] and |AAB[ Lemma 4.5] (see also ASS]) and we skip here the details for the sake of brevity, 
only stressing that the crucial ingredient for (4.53) is the fact that ||Vp TF ||oo < C. 

The second inequality is a consequence of the equation for Q, which for x£ Q, and y € Qj, i ^ j, 
can be rewritten 

V x p TF (z) ■ V x S(x, y) - p TF (x)A x g(x, y) = 0, 
together with the Gagliardo-Nirenberg inequality (see, e.g., [Nj Theorem 1]) 



|V X <5|I <C ||A X <5 



1/2 



1/2 



where the suprema are taken with respect to x G Q l and for almost every y G Qj. Combining this with 
the inequality above we get HVx^H^ < C || 1 | ^ si nce p TF is uniformly bounded from below inside the 
support of n+. Using (4.53) one then obtains (4.54). □ 



5 Energy Lower Bound and Convergence of the Vorticity Mea- 
sure 

In this section we complete the proof of Theorem |2.1| by proving the following: 

Proposition 5.1 (GP energy lower bound). 

If £1 = Hq \ loge|, with Hq > fli as e — > 0, then 

EGP > £GP +/ TF| log£ |2 f 1 + o nog^s\\\ (51) 



In the course of the proof we will gather the estimates leading to the conclusion of Theorem 2.2 



particular the rigorous version of (2.29). The strategy we employ makes use of several classical techniques 
which have already been used in similar contexts, see, e.g., jAABl IABMI [CRY, IM1, R . It originates 
from the GL theory of type-II superconductors (see |SS2| and references therein). 

5.1 Preliminary Steps 

We start as usual with the energy decoupling 

E GP = E GP + f dre g (u), (5.2) 
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where we have denoted 

u := (5.3) 

9 

and the reduced energy density is 

e fl («) := \g 2 \Vu\ 2 - fig 2 ^ • (iu, V«) + ^ (l - M 2 ) 2 , (5.4) 

with, as usual, 

(m, Vu) = — (hV«* — ii'Vu) . 

For a lower bound it is possible to neglect the contribution to the reduced energy in the region M. 2 \B(R+) 
with 

R+ :=i? TF +e 1 / 4 . (5.5) 
Indeed, the only possibly negative term can be bounded below by a simple trick 

-Q I drg 2 • (iu,Vu) > -\ [ dr g 2 \ Vu\ 2 - Q 2 [ dr 5 2 |w| 2 (5.6) 

Jr 2 \b(r + ) 4 Jm 2 \b(r + ) Jr 2 \b(r + ) 



and in K 2 \ B(R + ) we can use the exponential smallness (see Proposition A. 6 ) of $ GP = gu to bound the 
second term and obtain 

/ dre 9 (w)> / dr e g (u) - 0(e°°) (5.7) 

where positive but useless contributions have been dropped. 

Note for later use that the trick (5.6) used on the whole of ]R 2 together with the trivial bound 

£[u] < 0, 

which follows by taking g as a trial state for £ GP , imply 

j^dr |i 9 2 |V U | 2 + ^(l-H 2 ) 2 |<C|loge| 2 , (5.8) 

because of the normalization of g 2 |u| 2 = |\& | 2 . 
We now define 



R+ 

2/ 



F(r) :=-f2 / dttg 2 (t), (5.9) 

J r 

which satisfies llogelV^F = g 2 ttv^ in B(R+) and F(R + ) = 0. We can then integrate by parts the 
momentum term in (5.7| to obtain 

/ dre g (u)=[ dr \\g 2 \Vu\ 2 + | loge|Fcurl(iu, Vu) + ^ (l - |w| 2 ) 2 ) . (5.10) 

We will soon construct vortex balls to bound this energy from below by the energy of vortices of u. 
There is however one more reduction we need to perform before: as it is well-known, the vortex balls 
construction is feasible only in the region where the matter density is large enough. We thus split the 
domain again, introducing for some constant Cbuik > to be chosen later on (i?buik is the radius appearing 
in the statement of Theorem 2.1 ) 

^bulk : = R+ — Cbulk^ 1 , -Rcut : = R+ — ^^bulk^ (5-11) 
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Note that we have now four radii in our construction: 



-Rbuik < -Rcut < -R TF < R 



and in particular 



^cut < R 



■TF 



err 1 , 



(5.12) 

which is a consequence of (5.5) and (5.11). Using (A. 5) and the explicit form of p TF (in particular the 
fact that it vanishes linearly on dB(R lb )), this implies 



9 2 {r) > C\ log el" 1 , for any r <= B{R C 



(5.13) 



and for some C > 0. 

We now perform a smooth cut-off at i? C ut by introducing regular radial functions Xin and Xout satisfying 



Xit 



Xout = 1 in 



Xin(f) = 1 for r < i?buik, Xin(r) = for r > R cut . 



Using (5.11) and recalling that we consider the regime Q = 0(| logeQ we can impose 

|V X in| <C|loger\ |Vxout| ^Cllogel- 1 . 



(5.14) 



(5.15) 



The point of performing this cut-off is that we choose Cb u ik small enough in (5.11) in such a way that 

1 



|loge||F| < - A g\ in B(R+)\B(R huii ). 



(5.16) 

This is an easy consequence of the definition of F ( |5.9[ ): for any r € B(R+) \ B(Rbuik) one has (recall 
that g is decreasing as discussed in the Appendix) 



\F(r)\=n ( 



dttg 2 (t) 



< -n g 2 (r) (Rl 



<ic bulk J? + |log £ |-VW, (5-17) 



so that the above inequality holds true if one picks, e.g., Cbuik < \R+ 1 ■ Using (5.16) and the trivial 
bound (the two line computation is in [CPRY3, Lemma 3.4]) 



we see that 



This leads to 



Xout 



/ dr e g (u) > dr < 

Jb(r+) Jb{r + ) I 



Xir 



|curl(m, Vu)| < |Vu|' 



-5 2 |Vw| 2 + |loge|Fcurl(TO,Vu) I > 0. 



^g 2 \ Vm| 2 + Fcurl(iu, V«) 



2\2 



+ TXoutSlV^ + Mi - |<) 



(5.18) 



5.2 Evaluating the Individual Energy of Vortices 

We now state two classical results allowing to spot the vortices of the condensate and evaluate their self- 
energy. The methods (vortex balls construction and Jacobian estimates) are classical tools by now and 
we thus skip the proofs. The reader can refer to [SS21 Chapter 4 and 6] for a self-contained presentation 
of the techniques and to |IM1( Proposition 4.1], |AAB( Proposition 4.1], |CRY| Proposition 4.2 and 4.3] 
for applications in contexts similar to the present one. The original references are [JeJ IJS2L ISaj . 

The first result isolates the possible zeros and phase singularities, i.e., vortices, of u in small balls. It 
is a consequence of the bound (5.8) and uses (5.13), which makes the reduction to B(R cu t) mandatory. 
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Proposition 5.2 (Vortex ball construction). 

There is a certain Eq > such that, for s < Eq there exists a finite collection {Bj}j^j := {B(a.j, Qj)}j eJ 
of disjoint balls with centers &j and radii Qj such that 



1. {r g B(R cut ) : 

2. J2l>3 = 1 loge| 
is J 



(r)|-l| > |loge|- 5 } cU ieJ Ei, 



Setting dj := deg{u,dBj} , if Bj C B(R cu t), and dj = otherwise, we have the lower bounds 



jf drg 2 |Vu| 2 > nldjlg 2 ^) |loge| (l - C 



log | log g | 
|loge| 



(5.19) 



for some fixed constant C > 0. 



Note that the error term O (log | log e | ) in (5.19) has two origins: one is the use of the growth and 



merging method [SallJej . It could probably be improved by using more refined versions of the technique, 
as presented, e.g., in |SS3j . The other is the approximation of g 2 in each vortex ball by its value at the 
center using (A. 5) and the regularity of p TF . The second error is much smaller than the first one and is 



thus absorbed in the 0(log | loge|) term. It is also much more intrinsic to our setting since it originates 
from the inhomogeneity of the matter density. 

Remark that there is some freedom in the definition of the vortex balls: we could take any negative 
power of |loge| at point 1 of the statement above, obtaining ball families such that the sum of the 
radii would be bounded by an arbitrary negative power of | log e| . This would not affect the order of 
magnitude of the remainder in ( 5.19[) and would improve some other remainders that will appear in the 
sequel. Since the error we make in (5.191 will necessarily be present in the final result, we see no point 



in over-optimizing other estimates and thus stick to the concrete choice made above (which by the way 
is the same as in [IMTI Proposition 4.1]). 



We denote 



F[u] 



dr 



B(R+) 



,9 2 |V«| 2 



(1 



(5.20) 



It is well-known |JS2j that this quantity can be used to control the vorticity measure (recall our normal- 
ization choice) 

fi = |loge| _1 curl(zM, Vu) , (5.21) 
which is the content of the following result. 

Proposition 5.3 (Jacobian estimate). 

Let <f> be any piecewise-C 1 test function with compact support 



supp(0) C B{R hu 



ik) 



Let {Bj}j GJ := {B(slj, fj)}^ 7 be a collection of disjoint balls as in Proposition 
deg{u, dBj}, ifBj C yB(i? cu t), and dj = otherwise, one has 



5.2 



Setting dj 



2^ 2lrd ^ 



log el 



B(-Rbuik) 



<C|loge|- 9 ||V0|! 



L~(B(i?. bulk )) 



T\u\ 



(5.22) 



where the constant C is independent of <fi. 
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Taken together, these results imply 



f dr \ X in (lg 2 \Vu\ 2 + \loge\Fcurl(iu,Vu))} > i f dr Xin .g 2 |V«| 2 

JB(R + ) I \ Z / J 1 JB(R + )\U teI Bi 

+ ^2 7 r|logE| Xin (a i ) (IWiaA + diFfa)) (l - ) - OQ logs] log | loge|), (5.23) 



where the bound (5.15) has been used. The last remainder terms contains two contributions: the smaller 
one is due to (5.221 can be easily seen to be of order 0(\ log e | ~ 7 ) by using the a priori bound (5.8 1 on 
F[u\. The second one is due to the factor —Clog | log ej | loge| _1 multiplying the (negative) vortex energy 
gain £27r|loge|d 4 F(a,-), i.e., 

log | log £ | / dr Fcurl (iu, Vzt) = — Sl log | log er| / dr g 2 (iu, Vu) 
Jb(r+) Jb(r+) 



and it can be estimated precisely as in (5.6), again with the help of (5.8), yielding an error term of order 
|loge|log|loge|. 

The purpose of the next section is to prove a lower bound to the first term of the right-hand side of 
this inequality. As we will see, it is this term that contains the energy due to the interaction between 
vortices. 



5.3 Evaluating the Interaction Energy of Vortices 



We now estimate the first term of the right-hand side of (5.23). Following [ABMj |Rj we first regularize 
the vorticity measure /x, by introducing the modified current 



(*«,V«), in B(R cut ) \ (J 
0, otherwise. 



(5.24) 



The regularized vorticity measure jl is the one naturally associated with j (rcnormalized as in (2.13)): 

ft := |loge| -1 currj. (5.25) 
It is more regular than the original /i in the sense that we have removed the singular part of the phase 



generated by vortices by neglecting in (5.24) the part of the current inside the vortex balls. There is still 



of course a line singularity along the boundary of the balls because the current j goes brutally to there, 
but this is much less problematic that the vortices lying inside the balls, which are point singularities. 

The following lemma shows that the regularized vorticity measure is close to the original one in the 
(C*)* norm, allowing to control the error term due to the replacement of [i by its regularization \x. 

Lemma 5.1 (Approximation by a regularized vorticity measure). 

There exists a constant C > such that, for any piecewise-C 1 test function 4> with compact support 



supp(0) C B(Rcut), 



we have 



(/i -fi)<f> 



B(-Rout) 



<C|loger 8 ||V<; 



li-(B(Hout)) 



(5.26) 
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Proof. By the definitions (|5.21| and (|5.25|) 

/ (m-m) 



|loge| 



JB(R oni ,) 



Then, using (5.13), 

drV-^-j 



< II Vc 



lL°°(B(fl cut )) 



jloge| 



dr V-'-f 



dr |u||Vu| 



< C|loge|||V0|| t)) / dr 5 2 | u ||V U | 

< C|bg E |||V0|| £ « (B(JJeut)) (<5 / drg^ + S- 1 [ dr g 2 \Vu\ 2 ) 



for some 5 > 0. We now use the upper bound g 2 |it| 2 = l^ 013 ! 2 < C (see Proposition A. 6) and the upper 

to obtain 



bound on . £>j of Proposition 



5.2 



/ dr.g 2 | U | 2 < C|loge|- 

JUjsjBj 



20 



Recalling the bound (5.8), we also have 

/ dro 2 |Vu| 2 < C|loge| 2 
Js. 2 

and, choosing <5 = | log el 11 , we conclude the proof. 
We now define hjx as the unique solution to 

f-V (^V^ A ) = ft, in B(iJbuik), 
\/i A = 0, on dB(i? bulk ), 

and claim that the following holds: 

Lemma 5.2 (Vortex interaction energy). 

Let {Bj}j eJ := {S(aj, j &e a collection of disjoint balls as in Proposition 



□ 



(5.27) 



/ dr X i„.g 2 |Vu| 2 > (l-C|log £ r 5 )|loge| 2 / dr^V/ip 

Je(i?. + )\u 3£J e J JB(flbuik) 3 



(5.28) 



/or some given constant C > 0. 

Proof. The proof is the same as in [R] Lemma 3.3], only a bit simpler because we work in a simply 
connected geometry. We recall that \m = 1 in S(-Rbulk) \ UjGjBj an d |u| is close to 1 there, according to 
point 1 in Proposition 5.2 Thus 



B(fl+)\U 3 



dr Xi „g 2 |V M | 2 > f dr 5 2 |Vu| 2 

ijBj JB(ii bulk )\U je jB 3 



> (l-C|loge|- 5 ) f dr o 2 |u| 2 |Vu| 2 

"'B(i?. bulk )\U je . 7 B J - 

> (l-Cllogep 5 ) f dr.g 2 |(m,Vu)| 2 



(l-C|log £ r 5 ) f dro 2 |J 

JB(R^u) 



2 |~|2 

ui g 

B(Rhx 
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curl | logep 1 ] + -^V 1 ^ = 0, in H' 1 (6(i? bulk )) 



By definition 



so that there exists / € H 1 (B(Rb u ik j) satisfying 

j = -|log£|iv x ft A + V/. 
9 

Using the fact that hp is constant on dB(Rhuik), we have 

/ drV ± /i A -V/ = 

by Stokes' formula, and thus 

/ dr ff 2 |j| 2 = |log £ | 2 / d r4|V/* A | 2 + / dr ff 2 |V/| 2 

and there only remains to drop the last positive term to conclude the proof. 

5.4 Completion of the Lower Bound Proof 

Injecting ( |5.28| ) in ( |5.23 1, we have 



□ 



/ dr( X inf^.g 2 |V U | 2 + |log £ |Fcurl(z W ,V U )H> (l-C|log £ |- 5 )|loge| 2 f dr^|Vfc A 

JB(R+) L \ Z /J ^B(flbulk) z 9 

+ M logelXinK) Q|d,|.9 2 («,) + diFiaA (l C ^"^ ) - OQ loge\ log | loge|) 



and because of Proposition A. 4 we can replace g 2 and F by p TF and F TF respectively and absorb the 
new remainder in the old ones. Using in addition Proposition |5 . 3| we can also replace the sum over vortex 
balls and obtain 

/ dr(xinf^.g 2 |V U | 2 + |log £ |Fcurl(z W ,V U )U>(l-C|log £ |- 5 )|loge| 2 f dr-^\Vhtf 
JB(R+) V \ l J) JB(R bu m) Z P 

, log I log £ 



!loge| 2 1-C 



|loge| 



B(R C 



Xin 1 5 P ImI+^V) -0(|loge|log|loge|). (5.29) 



We can apply Lemma 5.1 with either <fi = XinP TF or <j> = %i n F TF , to deduce 

/ dr Xi „f^ 2 |V U | 2 + |log £ |Fcurl(z M ,V W )) > (l - C| log^ 5 ) | loge| 2 f dr |Vfy 

/ JB(R bum ) Z P 

,log|loge| 



|loge| 2 1-C 



| log 



p iJ, |£|+F"/i -0(|loge|log|log £ |). (5.30) 



Now the definitions (5.9) and (5.11) guarantees that inside B(R cut ) \ B(Rhu\k) 
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by the same computation as in (5.17), and we can thus drop one more positive term to arrive at (recall 
that Xin = 1 in S(-Rbulk)) 



/ drxi„ ( x.g 2 |Vu| 2 + |loge|Fcurl( iM ,Vu) ) > 
Jb(r+) \ z J 



loge|^ l-C 



,log|loge| 



! loge| 



B(-Rbuik) 



2p 



' |V^| 2 + ^ TF |A|+F TF A) -0(| loge| log | log £ |). (5.31) 



,TF I 



Gathering (5.2), (5.7|, (5.18) and the above inequality, we thus conclude 

, log | log e 



E GP >E 



GP 



loge| 



l-C- 



| loge| 



B(flbuik) 



2p 



Wl V %| 2 + ^ TF lAl+^ TF A 



0(| loge| log | loge|) (5.32) 



and one can recognize the renormalized energy I TF of the regularized vorticity measure jj, on the right- 
hand side, except that the domain is £>(i?buik) instead of B(R TF ). However it is clear from Theorem 



3.1 that B(Ri,) C S(i?buik) for any fixed fio and e s mall enough since i?* does not depend on e and 
-Rbuik — > R TF when e —> 0. We can thus invoke Remark 3j2 to deduce that the infimum of the renormalized 
energy appearing on the right-hand side of (5.32) (taken with respect to jl) is precisely |loge| 2 / TF , so 
that one has 

E GP > ^GP + /TF | log£ |2 _ q ( | log£ | bg | bg£ | ) ) 



which concludes the proof of Proposition |5.1| and thus of Theorem |2.1 
5.5 Convergence of the Vorticity Measure 



It is now a short way to the proof of Theorem 2.2 We apply the stability estimate ( 3.8 1 to the renormalized 
energy appearing in (5.32), with integration domain rest ricte d to S(i?b u ik)- Since such an energy has the 
same infimum and minimizer as I TF (see again Remark 3[2 ) , 



E GP > £GP 



| loge| 



r TF 



9„TfI V ^A-M* 



l-C 



log | log e | 
| loge| 



-0(e°°) 



where for any measure v we denote now by h v the solution of 



\ w Vh v )=v 1 in B(R hulk ) 



[h v = 0, ondB(R hu i k ). 
Combining this with the energy upper bound of Proposition |4.1[ we obtain 



^logjloge, ; 



„TF 1 



(5.33) 



(5.34) 



which clearly implies that 
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Indeed the estimate (5.34) together with (5.33) imply that, for any differentiable function <f> with compact 
support strictly contained in B(R cu t), 



e(fl bu i k ) 



B(flbuik) P 



< c 



log I log £ I 

I logel 1 / 2 



1/2 



B(flbuik) P 



TF 



1/2 



(5.36) 



by the Cauchy-Schwarz inequality. Note that we have also vindicated the claim contained in Remark |2|3| 
that the regularized vorticity measure can be estimated in a better norm that p. 



Finally, combining (5.35) with Lemma 5.1 we obtain 

(p - p„)<j) 



sup 

0GCi(B(i?. bu ik)) 



B(fl bu „ 



B(i? bu ik) P 



•TF 1 



1/2 



< c 



log I log £ I 
| log e| 1/2 



1/2 



li°o(B(JibuJv) 



where /i is defined in (2.13). 



Appendix 

In this Appendix we collect some useful but technical estimates. Most of them have already been proven 
in very similar contexts so, for the sake of brevity, we keep the proofs as short as possible, referring to 
other papers for full details. 

We start by investigating the properties of g, i.e., the minimizer of £ GP . By standard arguments one 
can show that g is unique, smooth, positive, radial and decreasing and it solves the variational equation 

- a" - r- l g' + r s g + 2 £ - V = e- 2 X GP g, (A.l) 

for rel+. 

The evaluation of £ GF on a suitable regularization of \J p TF as well as the trivial lower bound E GP > 
E TF yield the estimate 

£ GP =£ TF + 0(|log £ |), (A.2) 

which in turn implies the convergence of g 2 to p TF and a similar estimate of the difference between the 
chemical potentials: 

||S 2 - P TF || L2(R2) < 0(£^og^), A GP - A TF < O(e^\logi\). (A.3) 



A very simple investigation of the variational equation (A.l) shows that the maximum of g is bounded 
by the chemical potential A GP and thus by a constant, i.e., 

Ilffllz-^) < O(l). (A.4) 
The next result is an improvement of the estimate of the difference g 2 — p TF : 



Proposition A.4 (Pointwise estimate of g 2 ) 

Let = flo \ loge| with J7q > fixed, then 



sup 

r6[0,i? TF -e 2 / 3 | log £ | 2 / 3 ] 



\g 2 (r)- p TF (r)\ <0(e 2 / 3 |loge| 2/3 ). (A.5) 
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Proof. The result has already been proven in a slightly different setting in [CRY,, Proposition 2.6] (sec 
also [AABI UMlj for similar results) but we repeat the crucial steps for convenience of the reader. 

The result is obtained through a local analysis of the variational equation ( A.l ), which can be rewritten 
in the form —Ag = 2e~ 2 (p — g 2 )g, with 



p(r) := 



1 



A GP - 



(A.6) 



Note that in the ball B(R TF — ^e 2 / 3 1 log e| 2 / 3 ) one has the lower bound p(r) > Ce 2 ^ 3 \ log e| 2 / 3 , since 
\p(r)- P TF (r)\ <0(sV\^i\) « p^ F (r) ^O(s^\\0ge\ 2 /% if r < R TF - \e 2 / 3 \ loge| 2 / 3 , (A.7) 



by \k^. 

Now we consider an interval [ro — <5, ro + 5] with 

S= l £ 2/3| log£ |2/ 3) 

for some |e 2 / 3 | loge| 2 / 3 < rr, < R TF — e 2 / 3 | loge| 2 / 3 . In this region it is possible to construct explicit 
super- (see [CRTl Eq. (2.36)] or [Ml Proof of Proposition 2.1]) and sub-solutions (see [CRY] Eq. (2.47)] 
or |Serj ) to (A.l), which provide the estimates 

p (r + S) (1 + 0(e~)) < 5 2 (r ) < p (r - 6) (1 + 0( £ °°)) , 

where we have also exploited the fact that p is decreasing. It remains now to replace p(r ±d) with /5(r ), 
which yields an error of order £ 2 / 3 ] loge| 2 / 3 and use (A.7). □ 

Another important property of g is its exponential decay outside the support of p TF , i.e., for r > R TF : 
Proposition A. 5 (Exponential decay of g). 

Let Q be as in Proposition \A.4\ then there exists two constants c > and C < oo such that for any 

r >i ? TF +£ 2/3 



g 2 {r) < C£ 2/3 |loge| 2/3 exp 



c(r - R rF ) 

£ 2/3 



(A. 



Proof. It suffices to notice that W(r) := cie 1 / 3 ! logel 1 / 3 exp{— C2(r — i? TF )/e 2 / 3 } is a supersolution to 
(A.l ) for r > R TF + e 2 / 3 and for some constants c%, C2 with C2 > and c\ < oo: Indeed at the boundary 
by (A.5) 

W(R TF +e 2 / 3 ) = c 1 £ 1 / 3 |loge| 1 / 3 e- C2 > g(r), for r > i? TF , 

if we pick ci — 0{\) large enough. Moreover W satisfies — AW + c 2 e~ 4 ^ 3 W > 0. On the other hand 
(A.l) yields — Ag + Ce~ 4 ^ 3 g < because p < —Ce 2 ^ 3 in the same region, so that IF is a supersolution 

□ 



for C2 equal to v C > and the result is proven. 

We now turn to properties of ^P GP , that we summarize in 
Proposition A. 6 (Pointwise estimates of , 3> GP ). 

Let f2 be as in Proposition \A.J\ then there exists two constants c > and C < oo such that for any 
r>R TF +e 2 / 3 

|* GP (r)| 2 <Ce 2 / 3 |loge| 2 / 3 exp/ c(r " ^ 



r2/3 



(A.9) 



Moreover, there is a constant C such that |$ (r)| < C uniformly in M 2 . 

Proof. The proof is essentially the same as that of Proposition |A.5| plus some standard trick to extract a 
differential inequality for |\1/ GP |. The reader can consult |AAB1 Proposition 2.5], |IM11 Proposition 3.2] 
or |CRY| Section 2.2] where similar results are proved. □ 
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